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Abstract 

Background: Lutzomyia gomezi (Nitzulescu, 1931) is one of the main Leishmania {Vianna) panamensis vectors in 
Panama, and despite its medical significance, there are no population genetic studies regarding this species. In this 
study, we used the sequences of the mitochondrial gene cytochrome b/start of NADH1 and the nuclear elongation 
gene a-1 in order to analyze genetic variation and phylogeographic structure of the Lu. gomezi populations. 

Methods: A total of 86 Lu. gomezi individuals were captured in 38 locations where cutaneous leishmaniasis 
occurred. DNA was extracted with phenol/chloroform methods and amplification of genes was performed using 
PCR primers for mitochondrial and nuclear markers. 

Results: We found a total of 37 and 26 haplotypes of mitochondrial and nuclear genes, high haplotype diversity (h) 
for all three populations were detected with both molecular markers. Nucleotide diversity (n) was estimated to be 
high for all three populations with the mitochondrial marker, which was opposite to the estimate with the nuclear 
marker. In the AMOVA Ost recorded moderate (mitochondrial) and small (nuclear) population structure with 
statistical significance among populations. The analysis of the fixation index (Fst) used to measure the 
differentiation of populations showed that with the exception of the population located in the region of Bocas del 
Toro, the other populations presented with minor genetic differentiation. The median-Joining network of the 
mitochondrial marker reveled three clusters and recorded four haplotypes exclusively of localities sampled from 
Western Panama, demonstrating strong divergence. We found demographic population expansion with Fu's Fs 
neutrality test. In the analysis mismatch distribution was observed as a bimodal curve. 

Conclusion: Lu. gomezi is a species with higher genetic pool or variability and mild population structure, due to 
possible capacity migration and local adaptation to environmental changes or colonization potential. Thus, 
knowledge of the genetic population and evolutionary history is useful to understand the implications of different 
population genetic structures for cutaneous leishmaniasis epidemiology. 
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Background 

Cutaneous leishmaniasis is the most common form of 
leishmaniasis reported in the Republic of Panama; its 
clinical manifestations range from minor lesions to severe 
skin ulcers [1-3]. In Panama it was first recorded in 1910 
[4], an overall of 15 cases occurred during 1910-1944; 
but there has been a sharp increase since 2000 [5]. 
According to the Epidemiology Department of the Health 
Ministry of Panama, a total of 26,163 cases of cutaneous 
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leishmaniasis occurred during 2000-2010 and the distribu- 
tion of occurrence of CL was in regions of the provinces of 
Bocas de Toro (28%), Panama West (20%), Code (17%), 
Colon (11%), Panama East (5%), Veraguas (3%). The disease 
is mainly caused by Leishmania {Vianna) panamensis and 
Lutzomyia ylephiletor, Lu. sanguinaria, Lu. panamensis, 
Lu. trapidoi, Lu. gomezi were identified as the transmission 
vectors [6-8]. 

In Panama, recent research has show that Lutzomyia 
{Lutzomyia) gomezi (Nitzulescu, 1931) is the most abun- 
dant species with wide geographical distribution and 
their abundance has been associated with cases of clinical 
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CL acquired in households of rural communities [9]. 
Alongside Lu. gomezi other Lutzomyia species vectors 
have been caught in the same place and associated with a 
focus of cutaneous leishmaniasis infection in other 
American countries, incriminating this species with high 
potential in the transmission of the disease [2,10-14]. 

Lutzomyia (Lutzomia) gomezi in Central America is 
present from Mexico to Panama and Trinidad Tobago 
(Caribbean), also recently reported in Guatemala [14]. In 
South America, it has been reported in Colombia, 
Venezuela, Ecuador, Peru, French Guiana and Brazil, 
specifically in Amapa, Acre, Para, Mata Grosso, Goias, 
Bahia, Maranhao and Rondonia [15-17]. 

In general, its natural environment is a humid and 
dark place, such as nests, rock crevices, animal burrows 
and tree bark in the tropical rainforest [2]. However, it is 
also reported in forest gaps and the canopies were light 
availability produces humid change [2,18]. The foraging 
and blood feeding behavior occurs at twilight and at 
night from 18:00-20:00 hours [18]. Nevertheless, the 
deforestation and loss of natural habitat have caused its 
adaptation to peridomestic areas, perhaps changing the 
hourly activity and feeding on a large variety of domestic 
animals [19,20]. In the focus of leishmaniasis, this spe- 
cies is predominant near peridomestic or outdoor areas 
more than indoors households, which is difficult to con- 
trol [2,19,21-23]. 

Despite of the potential significance of Lu. gomezi as a 
Leishmania vector, few researches have targeted genetic 
aspects or interrelationship of host-vector species. For 
instance, [24] detected a natural infection of Lu. gomezi 
with the Le. braziliensis in Venezuela, in order to establish 
a methods for determine the circulation of Leishmania 
parasites in leishmaniasis endemics areas. On the other 
hand, [25] analyzed the changes to the primary and 
secondary structures of tRNA ser of the species; Lu. 
trinidadensis (Oswaldoi group), Lu. (Psychodopygus) 
panamensiSy Lu. (Micropygomyia) cayennensis cayennensis, 
Lu. dubitans (Migonei group), Lu. {Lutzomyia) gomezi, Lu. 
rangeliana (ungrouped) and Lu. evansi (Verrucarum 
group) for taxonomic purposes, considering that morpho- 
logical identification can be difficult. Also, [26] detected 
with ITS-1 a pool of Lu. gomezi infected with Le. naiffi in 
Panama, the first report for the country which prompted 
several hypotheses on the introduction of this parasite into 
this country. 

Attempts to understand the role that arthropod vectors 
play in disease dynamics and pathogen transmission of 
leishmaniasis, several studies over genetic population of 
sandfly vectors has been performance in the Latin America, 
Iran, Turkey, Palestine, Israel and Egypt [27-30] . Many of 
these studies are focused in two principal species, Lu. 
longipalpis and Phlebotomus papatasi to determinate 
their genetic variation, structures and differentiation of 



populations [27,28,31,32]. Thereby, molecular evidence for 
divergence of vectors has been found and an assessment of 
the impact on leishmaniasis epidemiology. For instance, 
the taxonomic status of the Lu. longipalpis complex has 
been fundamental to the understanding of leishmaniasis 
epidemiology. Symptoms observed for transmission of 
Leishmania infantum chagasi by Lu. longipalpis in Brazil 
and Colombia results in visceral infections, whereas the 
transmission of the same parasite by Lu. longipalpis in 
Costa Rica results in non-ulcerative lesions [33,34]. 

Thus, it becomes important to analyze the population 
genetics and demographic history of this vector species 
to assess the genetic pool, gene flow and colonization po- 
tential, consequently known as the population changes 
and genetic variation affecting the vector competence 
and resistance to insecticides [35,36]. On the other hand, 
the population genetics analysis can provide information 
about migration capacity, reproduction, and adaptation 
to the conditions of the new habitat, favoring the emer- 
gence of vector diseases [37]. In addition to this, the 
determination of cryptic vector species and their ability to 
transmit pathogens are even more relevant to understand 
their implications in the epidemiology of vector diseases 
and to suggest appropriate and effective prevention and 
control programs without any environmental risks [36]. 

Due to the lack of population genetics information of Lu. 
gomezi species in Panama and the Americas as a whole, 
the goal of this research was to evaluate and compare the 
intra- and inter-population variability of Lu. gomezi from 
localities with a high incidence of cutaneous leishmaniasis 
in Panama. Also to identify the barriers that may influence 
gene flow among Lu. gomezi populations. Moreover, to 
infer the historical processes that defines actual geographic 
distributions of the species. These results may contribute 
on the knowledge of leishmaniasis epidemiology and im- 
prove the development of focal or large-scale programs 
control for leishmaniasis vectors in Panama. 

Methods 

Study sites 

The study was performed in the Republic of Panama be- 
tween the coordinates 7°11'-9°39'N and 77°10'-83°03' 
W. The most prevalent climatic regime in Panama is 
tropical humid, with the dry season (January-March) pre- 
senting an average temperature of 31.5°C and relative 
humidity of 75%. In other months (April-December), the 
average temperature is 27°C and the relative humidity 
averages 90%. The vegetation in Panama varies according 
to climate zones and consists mainly of tropical humid 
forest or savannahs resulting from agricultural activity. 

Sandfly collections and identification 

Lu. gomezi samples were collected from the thirty-eight 
localities in Panama, during the dry season (January- April) 
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and rainy season (May-June) of 2010. Thirty-four sites 
were chosen, representing incidence with cases of cutane- 
ous leishmaniasis in humans during the 2006-2009 period, 
according to data from the Department of Epidemiology of 
the Ministry of Health of Panama; four sampling sites were 
in Panamanian tropical forests. The geographic distribution 
of the samples locations are indicated in Figure 1. 

The characteristic of the thirty-four collecting sites 
according to the landscape was fragmented environment 
(secondary vegetation and areas used for agricultural 
activities) or rural (an area without a large concentration 
of people). The Coiba National Park, San Lorenzo Park, 

011 Pipeline Road Park, and the Serrania Filo del Tallo 
Park represent natural reserves (forests) protected by the 
National Environment Authority of Panama (ANAM). 
The samples localities are listed in Table 1. 

Specimens of Lu. gomezi used in the study were 
caught using CDC light traps and octenol solution uti- 
lized to attract hematophagous insects [38]. The traps 
were positioned in the fragmented and rural locations, 
near peridomestic areas, with the presence of domestic 
animal (dogs, chickens, pigs, and cats) and ornamental 
vegetation. In the forest, traps were placed next of stone 
crevices and tree barks. 

A total of nine CDC light traps were exposed at a 
height of 1.5 m in each locality and were installed at ap- 
proximately 50 m intervals. The collection period was of 

12 h (18:00 h-06:00 h) for two consecutive days at each 
point. The specimens were sacrificed with chloroform 



and stored in 95% ethanol at -20°C. Adult Lu. gomezi 
were identified by morphology and genital structure 
using the identification keys of [16]. 

DNA extraction and amplification 

The DNA of 86 individual sandflies was isolated by a 
standard phenol/chloroform technique, precipitated 
with ethanol and resuspended in Tris-EDTA (TE) buffer 
according to [39], with some changes and stored at -20°C 
until used. 

Two mitochondrial and nuclear genes were analyzed 
to estimate the structure and geographic genetics 
within and among populations of Lu. gomezi. The 
region of mitochondrial genes analyzed was final Cytb 
gene, inter-genetic region IGS-1, tRNA-Ser, inter-genetic 
region IGS-2, and start of the NAD HI (CB3-NIN), and 
amplified using primers CB3-PDR: 5 ' C AYATTC AAC 
CWGAATGATA37N1N-PDR:5'GGCAYWTTGCCTC 
GAWTTCGWTATGA3 ' [27]. Nuclear region elongation 
factor alpha- 1 (EF a-1) coding an intron region were amp- 
lified with primers EF-F03 5 CCTGGACATCGTGATTT 
CAT3 7EF-R04 5 ' AGTGCTTCGTGGTGTAT(C/T)TC3 ' 
[40]. 

PCR amplification of genes were performed in a volume 
of 25 ul containing lx buffer, 1.5 mM MgCl2 ++ , 0.2 mM 
of each dNTPs, 0.4 uM of each primer, 1.5 units of Taq 
polymerase and approximately 2 ul of DNA («50 ng/ul) 
template. The details of thermal cycling conditions were 
outlined by [27] for CB3-N1N and [41] for EF a-1. A 




Figure 1 Map of the Panama isthmus and distribution of the sampling sites of Lutzomyia gomezi. The arrow indicates the Cordillera 
Central and line (-) represents the Panama Canal. The numbers inside the square symbolize four natural reserves (forest) protected by the 
National Environmental Authority (ANAM), used as control. 
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Table 1 List of sampling sites with information on landscape features, CB3-N1N and EF a-1 haplotypes identified in 
each location 



Pnmiiatinnc 

1 Vs k/ Wl 1 G L 1 13 


1 orations 


Coordinates 


Macro 
Habitat 


Micro 
Habitat 


CB3-NIN 
haplotype 


EF a-1 
haplotype 




1 -Nance Valle Risco 


9.255 


-82.470 


Fr 


P 


H1-H2 


H1-H2 


Western populations 


2- Bisira 

3- Parque Nacional Coiba 


8.899 
7.627 


-81.862 
-81.730 


R 

Fo 


P 
S 


H3-H4 
H6 


H1 

H4 




4-Altos de Piedra 


8.515 


-81.087 


Fr 


P 


H5 


H1-H3 




5-Villa del Carmen 


8.800 


—80.552 


R 


P 


H7 


H1 




6-Coclesito 


8.811 


-80.550 


R 


P 


H9-H11 


H1 




7-Molejon 


8.767 


-80.51 3 


Fr 


P 


H8-H9-H10 


H1-H12 




8-Cutevilla 


8.772 


—80.487 


Fr 


P 


H10-H1 1-H12 


H1-H24 




9-Vaquilla 


8.699 


-80.1 95 


Fr 


P 


H10-H12 


H1-H23 




10-Chirigui Arriba 


8.663 


-80.1 87 


R 


P 


H13-H14 


H3-H22 




1 1-Quebrada Leona 


9.1 77 


-80.1 39 


Fr 


P 


H15-H16-H17 


H1-H17-H18-H19 




12-Cuipo 


9.090 


-80.051 


R 


P 


H10-H15-H22-H23 


H1-H6-H9-H20-H21 




13-Providencia 


9.207 


—79.999 


R 


P 


H9-H15-H21 


H1-H9 




14-Achiote 


9.226 


—80.030 


R 


P 


H12-H19-H20 


H1 


v — tr l l L l a l |JU|JUlcl UUl lo 


15- Pina 

16- Union Piha 


9.244 
9.273 


-80.040 
-80.021 


Fr 
Fr 


P 
P 


H10-H18 
H10-H15-H22-H23 


H1 
H1 




17-Parque Nacional San Lorenzo 


9.351 


-79.973 


Fo 


s 


H10 


H1 




18-Parque Nacional Camino del Oleoducto 


9.1 19 


-79.700 


Fo 


s 


H15 


H1 




19-Altos de Divisa 


9.1 19 


-79.693 


R 


p 


H10 


H1 




20-Cacao 


8.748 


-80.017 


Fr 


p 


H10-H28 


H1-H5 




21 -Trinidad 


8.808 


-80.019 


Fr 


p 


H29 


H1 




22-Vista Alegre 


8.808 


-80.014 


Fr 


p 


H15-H23 


H1-H11 




23-Valdeza 


8.789 


-79.965 


R 


p 


H10-H27 


H1-H9-H10 




24-Caimito 


8.739 


-79.937 


R 


p 


H23-H29 


H6 




25-Limon 


8.697 


-79.904 


Fr 


p 


H24-H25-H26 


H1-H7-H8 




26-Ollas Arriba 


8.804 


-79.912 


R 


p 


H10 


H1-H9-H12 




27-Madroho 


9.284 


—79.1 34 


Fr 


p 


H10 


H14 




28-Gato Real 


9.267 


-79.1 1 9 


Fr 


p 


H15-H18-H30 


H1-H14-H15-H16 




29-Buenos Aires (Chepo) 




— / O.o I / 


Fr 


p 


H10-H15 


H1 




30-Torti 


8.981 


-78.571 


Fr 


p 


H10-H15-H31 


H1-H13-H14 




31-Arimay 


8.694 


-78.146 


R 


p 


H15 


H1 


Eastern populations 


32- La Cantera 

33- Nicanor 


8.640 
8.544 


-78.169 
-78.035 


R 
Fr 


p 
p 


H15-H32 
H10-H33-H34 


H1-H13-H14 
H1-H13 




34-Buenos Aires (Darien) 


8.530 


-77.961 


R 


p 


H10-H15 


H1-H8-H9 




35-Parque Serrania Filo del Tallo 


8.465 


-77.993 


Fo 


s 


H10 


H1 




36-Bijagual 


8.459 


-78.012 


Fr 


p 


H15-H23-H35 


H1-H13 




37-Rio Iglesia 


8.401 


-78.007 


R 


p 


H15-H36-H37 


H8-H14-H25-H26 




38-Puerto Obaldia 


8.669 


-77.429 


Fr 


p 


H15 


H1 



Fr: fragmented forest; R: rural; Fo: forest; P: peridomestic; S: selvatic. 



negative control was included in each PCR reaction and Sequencing and alignment 

visualization of PCR product using 1.5% agarose gel with a PCR products were purified and sequenced directly by 
100 bp ladder. the Macrogen INC., Seoul, Korea Sequencing Service in 
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both directions with the same primers used in the amplifi- 
cations. The sequences were compared to the nucleotide 
BLAST (Basic Local Alignment Search Tool) tool available 
in GenBank NCBI (National Center for Biotechnology 
Information) in order to verify the similarities between the 
sequences of Lu. gomezi stored in the database. The 
CodonCode Aligner software was used to edit the nucleo- 
tide sequences and identify the heterozygous as double 
peaks on the chromatograms of the nuclear region. The 
sequences were aligned with the Muscle option included 
in the MEGA 5.05 software [42]. The polymorphic sites, 
single substitution (singleton) and average base frequen- 
cies were computed with MEGA 5.05 software. 

Data analysis 

To carry out an analysis of population genetics (genetic 
diversity, structure and differentiation) to specimens 
collected in different localities, we clustered and defined 
samples in western, central and eastern populations from 
Panama (Table 1). Thus, the specimens collected in Bocas 
del Toro, Veraguas and Coiba represented the western 
populations, those collected from Code, Colon, and west 
of the Panama canal were central populations, whereas 
specimens from east of the Panama canal, Darien and San 
Bias corresponding to the eastern population (Figure 1). 

Genetic diversity 

We emphasized that heterozygotes detected in the EF a-1 
were treated with recombination model in PHASE option, 
incorporated in the DnaSP v.5. To estimate the genetic 
diversity we used the following indexes: haplotypes num- 
ber, haplotype diversity (h), and nucleotide diversity (tt) 
for each population. DnaSP v.5 [43] and Arlequin 3.11 
[44] software was used for this purpose. 

Genetic structure and differentiation 

The population structure was determined performing an 
Analysis of Molecular Variance (AMOVA) [45,46]. The 
Arlequin 3.11 [44] software analyzed the variance compo- 
nents and significance levels p = 0.01 using non-parametric 
permutations of 1,000 times. To estimate levels of genetic 
differentiation among the populations a pairwise compari- 
son test was performed. Non-parametric permutations of 
1,000 times and a significance p = 0.01 was used to esti- 
mate Fst-value [45] in Arlequin 3.11 [44] software. 

Network haplotypes 

A haplotype network for both genes was inferred 
through the median- joining network method [47], using 
the Network v.4.5 software. In order to obtain the values 
of posteriori probabilities among the clades, a Bayesian 
analysis was calculated using MrBayes 3.1 [48] program. 
MrModelTest [49] software was used to choose the best 
model of nucleotide substitution according to the AIC 



(Akaike Information Criterion). The model obtained by 
MrModeltest for CB3-N1N Bayesian inference analysis 
was HKY + I, with the proportion of invariable places (I) = 
0.7594. For EF a-1, the model obtained was K80 + 1, with a 
proportion of invariable places (I) = 0.6725. 10 million gen- 
erations of Markov Chain Monte Carlo (MCMC) and 
burn-in of 5 million were used. 

Demographic history of population 

In addition to this, Tajimas D [50] was used for estimate 
deviations from neutrality mutation or selective neutral- 
ity of these populations, P-value 0.01 was considered 
statistically significant. Also Fus Fs [51] were calculated 
for neutrality of mutation against population growth, 
hitchhiking and background selection, P-value 0.02 as 
significant. Both values Tajimas D and Fus Fs were 
calculated using the software Arlequin 3.11 [44] with 
permutations of 1,000 times. 

The mismatch distribution was calculated for CB3-N1N 
and EF a-1 as supplementary test to evidence of demo- 
graphic expansion of the population, through sudden 
expansion model [52]. In Arlequin 3.11 [44] we computed 
the sum of square deviations (SSD) and Harpending's 
raggedness index to test the goodness-of-fit of the model 
to the data. The time since population expansion was 
inferred using CB3-NIN region by means of the equation 
t = T 12 u, where tau (t) was obtained on the mismatch 
distribution outputs and a mutation rate (u) of 1.1-1.2% 
per million years in Heliconius erato [53]. 

Results 

Sequence characterization 

Overall, 86 sequences isolated of CB3-NIN and EF a-1 of 
Lu. gomezi individuals from 38 selected localities from 
Panama Isthmus were analyzed. Of these, an alignment of 
501 nucleotides was obtained of CB3-N1N with thirty- 
nine polymorphic sites and nineteen single substitutions 
(singleton). The average base frequencies were A = 37%, 
C = 15%, G = 10%, and T = 38% recording a sequence rich 
in A-T (75%). A fragmented size of EF a-1 was 597 bp, in 
this case 166 sequences were obtained by clone generating 
with the recombination model algorithm of DnaSP v.5 
considering the heterozygous individuals. Of those 
sequences, thirty-three polymorphic sites, six singletons 
and a larger proportion of G-C (53%): 23%T, 27%C, 23%A, 
27%G, were detected. 

Genetic diversity 

For the analysis of genetic diversity indexes both CB3-NIN 
and EF a-1 region demonstrated high haplotype diversity 
(h) for all three populations [Table 2]. While the nucleotide 
diversity (tt) for the mitochondrial region was estimated to 
be high for all three populations, nuclear nucleotide diver- 
sity was opposite to the estimate with mitochondrial, 
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Table 2 Genetic diversity data obtained with analysis of 



CB3-NIN and EF a-1of the Lu. gomezi 


CB3-NIN 


Western 


Central 


Eastern 




populations 


populations 


populations 


N° individual 


7* 


51 


28 


N° haplotype 


6 


23 


12 


h(d)±SD 


0.95 ±0.096 


0.89 ± 0.033 


0.80 ±0.061 


n(d) ± SD 


0.01 ±0.008 


0.01 ±0.005 


0.01 ±0.006 


EF a-1 


Western 


Central 


Eastern 




populations 


populations 


populations 


N° individual 


7* 


51 


28 


N° haplotype 


4 


18 


9 


h(d)±SD 


0.58 ±0.1 37 


0.58 ± 0.060 


0.63 ± 0.069 


n(d) ± SD 


0.003 ± 0.002 


0.004 ±0.003 


0.003 ± 0.002 



H(d): haplotype diversity; n(d): nucleotide diversity; SD: standard deviations* 
We highlight that the collecting sites of Bocas del Toro and Veraguas were 
difficult to access therefore have a much smaller sample size. 

showing a low value [Table 2]. The haplotype and the land- 
scape features per collection site are showed in [Table 1]. 

Genetic structure and differentiation 

In the AMOVA Ost estimated between the three popu- 
lations, CB3-N1N region recorded moderate and EF a-1 
small population structure with statistical significance 
among populations. The total genetic variation found 
within the populations was « 86% (CB3-NIN) and « 95% 
(EF a-1) as shown in Table 3. Pairwise comparison ana- 
lysis of the mitochondrial region revealed strong genetic 
differentiation with highly significance between Western 
Panama and Central Panama, whereas the Fst-value be- 
tween Western and Eastern Panama, showed a moderate 
and significant differentiation. However, when popula- 
tions were analyzed with nuclear region it revealed weak 
differentiation among populations of Western Panama 
and Eastern Panama. Also a moderate level of genetic 
differentiation was observed among the populations a 
Central Panama and Eastern Panama [Table 4]. 

Network haplotype 

The median- Joining network reconstructed with CB3-NIN 
region linked a total of 37 different haplotypes with three 
defined clusters (A, B, C). The lack of hierarchic structure 



is observed in the cluster "A", this clade consists of 28 hap- 
lotypes, the most frequent, H10, was found in eleven 
localities sampled in Central Panama and six in Eastern 
Panama. It also included Altos de Piedra and Coiba haplo- 
types from Western Panama. Cluster "B" is constituted by 
four haplotypes. Of these, H15 was collected from six sites 
in Central Panama and nine from Eastern Panama, with 
H23 being the most frequent [Figure 2A]. In cluster "C" we 
recorded four haplotypes exclusively of localities sampled 
from Western Panama, demonstrating strong divergence 
support by Bayesian analysis [complementary data]. Thus, 
the network is characterized by few mutational steps 
[Figure 2A]. 

On the other hand, analysis linked a total of 26 haplo- 
types of the EF a-1 region demonstrating a lack of genetic 
structure in the Lu. gomezi populations. The patterns of a 
star-like shape as similarly obtained with the mitochondrial 
region, showed haplotype HI as the most frequent in the 
center and several haplotypes surrounding it with one and 
three mutational steps [Figure 2B]. 

The Bayesian phylogenetic tree for the CB3-NIN region 
showed three clusters with a posteriori probability of 
0.99-1.0. Other minor clusters were supported by low a 
posteriori probability values. The cluster with a higher a 
posteriori probability corroborated the divergences obtained 
in the haplotype network, inferred by the median- joining 
method [Figure 2A and B]. The Bayesian inference of EF 
a-1 region did not show similitude either, this analysis per- 
formed for mitochondrial and nuclear genes are supple- 
mentary data attached as Additional files 1 and 2. 

We found 50% of the haplotypes in the peridomestic 
areas of fragmented environments, 40% in rural environ- 
ments, and 10% in forest environments. Although this 
information is additional and not the main focus of the 
research, the predominance of the haplotypes in these 
sites is important to highlight due to the anthropogenic 
impact on the landscape of Panama [Table 1]. 

Demographic history population 

The results of Tajima D test (CB3-NIN) were negative but 
not significant (-0.998; p > 0.01) that indicated neutrality 
deviation due to a demographic expansion or not neutral 
selection. Nevertheless, Fu's Fs values were significantly 
negative (-19.68; p < 0.02) detecting that a population 



Table 3 Analysis of molecular variance (AMOVA) based on CB3-NIN and EF a-1 sequence of the populations of 
Lu. gomezi 







CB3-NIN 






EF a-1 




Source of variation 


d.f 


(%) 
variations 


Fixation index 


d.f 


(%) 
variations 


Fixation index 


Among populations 


2 


13.7 


0 ST = 0.13* 


2 


4.5% 


0 ST = 0.04* 


Within populations 


83 


86.3 




163 


95.5% 





d.f = degrees of freedom; * P < 0.01 . 

All samples were grouped according to Table 1. 
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Table 4 Estimates of pairwise Fst of CB3-N1N and EF a-1 
between Lu. gomezi populations 



CB3-NIN 
EF a-1 



Western populations Central populations 



Western populations 






Central populations 


0.23 






0.02 




Eastern populations 


0.21 


0.10 




0.04 


0.12 



Significance (p < 0.05)*. 

expanded in the past. The assessment of demographic 
change with EF a-1 region using both Tajima D and Fu's 
Fs was significantly negative (-2.08; p < 0.01 and -13.65; 
p < 0.02) respectively. 

Considering that Fu's Fs are sensitive estimators for de- 
tecting demographic population expansion, we included 
the analysis of mismatch distribution for characterization of 
expansion detected with both genes. A bimodal curve was 
observed with CB3-NIN region and smoother curves with 
EF a-1 region. SSD value (0.011; p > 0.05) for CB3-NIN did 
not reject the sudden expansion model, while SSD value 
(0.049; p>0.05) for EF a-1 region was observed. The 
Harpending's raggedness index determined a good fit of 



data (Additional file 3). The time since expansion of 
the population was estimated to be approximately 
4,223,825 years ago corresponding to the Pliocene Epoch. 

Discussion 

The use of the mitochondrial and nuclear marker in ecol- 
ogy has provided enough information to allow contrast of 
the past population processes of species [54]. The high 
content of A-T (74%) of the mitochondrial gene analyzed 
in Lu. gomezi is similar to that detected in other sandfly 
species: Phlebotomus papatasi and Lu. evansi [30,55]. This 
suggests a high rate of substitutions as a result of natural 
selection, favorable or unfavorable to populations [56]. A 
higher proportion of G-C (53%) observed in the nuclear 
gene in Lu. gomezi was comparable to the proportion veri- 
fied in Anopheles gambiae (45%) and Culex pipiens (63%), 
characteristically of intron region [57,58]. 

The high haplotype diversity detected among popula- 
tions indicated a large polymorphism or gene pool of the 
Lu. gomezi. The existence of this polymorphism is im- 
portant for the survival and adaptability, considering that 
the natural habitat of Phlebotominae sandflies is 
destroyed because of anthropogenic factors [20]. How- 
ever, we observed discordance in the level of nucleotide 
diversity between CB3-NINI and EF a-1 sequences. This 




Figure 2 A haplotype network inferred by a median-joining method using CB3-NIN sequences (A) and EF a-1 sequences (B) of 
Lutzomyia gomezi. Circles represent different haplotypes and sizes that are proportional to haplotype frequencies. The length of the traces 
indicate number of mutations among haplotypes. Each color of circles showed the geographical populations: Western (red), Central (yellow), 
Eastern (blue). Values in parentheses showed the a posteriori probability of Bayesian inference. 
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may be because the mitochondrial gene is highly variable 
due to an elevated mutation rate and inheritence mater- 
nally [59]. While the EF a-1 is an intron region of low 
copy number of protein coding highly conserved [41]. 

The result obtained from AMOVA showed mild genetic 
structure but greater genetic variation within populations. 
The explanation is that gene flow occurs between the 
different populations, probably caused by the flying cap- 
acity of the species increasing their migratory potential. 
Similar results were found from the analysis of popula- 
tion genetics of Lu. longipalpis, Lu. intermedia and Lu. 
shannoni over a distance between 8 to 15 km, which sug- 
gested that individuals from sandfly populations were 
capable of flying more than 57 m, as previously recorded 
for sandfly species [2,60-64]. Thus, we hypothesized that 
Lu. gomezi is able to disperse over long distances, which 
could be gradually helped by the winds (adults) or passive 
movement of eggs and larvae. Studies across distant 
localities using other molecular markers are necessary to 
define vector capacity to disperse. 

Meanwhile, pairwise Fst-value demonstrated significant 
genetic differentiation with the CB3-NIN region of the pop- 
ulations from Western Panama and Central Panama. We 
think that this genetic differentiation observed is influenced 
by the central mountain range with altitudes of 2,468 to 
3,475 m. The mountain range represents a geographic 
barrier restricting altitudinal distribution of Lu. gomezi, 
estimated between 175-300 meters high [65]. Similarly, 
strong genetic differentiation was found between 
Anopheles albimanus populations from Costa Rica and 
Panama, also suggesting that the Central American 
mountain range is a physical barrier that limited the 
gene pool of malaria vector mosquitoes [66,67]. 

On the other side, slight genetic differences were de- 
tected from CB3-NIN and EF a-1 regions, among the 
sandfly populations from Eastern Panama and Central 
Panama (separated by the Panama Canal basin) could 
be explained because of the variety of climate and vege- 
tation in Panama. The climate of Panama is influenced 
by the Inter-tropical Convergence Zone, oscillating one 
year to other and from one region to another. This vari- 
ation allowed three types of tropical climate: Afi (very 
humid), Ami (humid), Awi (dry), according to Koppen 
Climate Classification [68]. The sampling locations from 
Veraguas, Code, Colon, and Panama (west of the province 
of Panama) registered precipitations and temperatures be- 
tween 3,001-3,500 mm and 22°C-24°C, respectively, while 
in the region of Eastern Panama, the amount of rain and 
the temperature reported varied between 1,000-1,800 mm 
and 26°C-27°C, respectively. 

As a result, the type of tropical climate of Panama pro- 
vides enough conditions to enable vegetation heterogeneity 
and refuge for a great diversity of species [69]. Lu. gomezi 
samples from Central Panama were collected from 



borderline locations of tropical forest vegetation of low 
lands. While the Eastern Panama populations were from 
areas of vegetation in regeneration with evident an- 
thropogenic impact. Thus, both climate and vegetation 
defined different ecotypes that can restrict gene flow due 
to environmental gradient. Climate changes were also 
used by [70] to justify the genetic differences found 
among the Brazilian populations of Lu. whitmani. 

The star-like shape of network haplotypes, both for CB3- 
NIN and EF a-1, reflected a rapid and recent expansion of 
Lu. gomezi populations. The haplotypes of CB3-NIN (H-10, 
H-15, H-23) and haplotype EF a-1 (H-l) are older alleles or 
ancestral haplotypes dispersed in the Isthmus. The network 
patterns in clusters A and B for the CB3-NIN region and 
network EF a-1 region could be explained due to the 
sandflies migration with their hosts, mainly mammals. 

Several mammal species (Canidae, Procyonidae, Sciuridae, 
Didelphidae, Dasypodidade) that currently are sandfly 
hosts established themselves in Panama during the north- 
south migration that took place during the closure of the 
Isthmus between the Miocene-Pliocene epochs [71-73]. 
We suggest this event as causal of the dispersal of Lu. 
gomezi because fossil records from the Miocene showed 
the presence of sandflies in fragments of mammal hair and 
in the microstructure of bird plumages in the Dominican 
Republic, suggesting their relationship with several hosts 
and possibly their transportation mechanism [74]. 

The H6 of CB3-NIN region from Coiba National Park 
(Veraguas) probably derived from the Hll detected in 
the populations of Central Panama, which indicates a 
possible migration of Lu. gomezi from the continent. 
Coiba was a continental island separated from South 
America in the Eocene, laying next to the Isthmus of 
Panama due to the decrease in the sea level in the last 
glaciation. This could have enabled the migration of wild 
animals and later, when the sea level increased, Coiba 
was isolated [75,76]. These factors could have isolated 
this haplotype, and it is also evidence that wild hosts 
could help disperse these vectors. 

The cluster B clearly showed the divergence in a west 
population (Bocas del Toro) and this strengthened 
significantly with the Fst-value [Table 4], we suppose 
that the population could have become geographically 
isolated. However other observations are necessary to 
establish phenotypic differences and corroborate this 
assumption. 

The fragmentation of the forests on Panama is 
reported since pre-Columbian times [77,78]. The 
prevalence of 10 (CB3-NIN) and 1 (EF a-1) ancestral 
haplotypes all types of environment and the percent of 
haplotypes found in the fragment and rural environ- 
ments could indicate the displacement of Lu. gomezi 
because of disturbed areas. However, it is necessary to 
establish the consequences of fragmentation forest in 
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the population genetics by increasing sampling in local 
forests. 

The neutrality tests Fu's Fs results were significantly 
negative, and mismatch distribution statistic and network 
haplotype support the assumption of population expansion 
in the past. Moreover, high haplotype diversity but low nu- 
cleotide diversity in EF a- region also suggests population 
expansion. 

Is it the fact that the closure of the Isthmus of Panama 
affected by climate change in the world was the main fac- 
tor that enabled the establishment of several species in 
America, allowing for faunal exchanges between North 
and South [79] . An alternative explanation of the popula- 
tion expansion is that conditions later to closure of the 
isthmus may have been favorable to its adaptation to the 
suitable habitat and different ecological conditions in 
Panama territory. According to our estimated date the 
Pliocene climate was much cooler, allowing the expan- 
sion of plant and animal species into new habitats [80]. 
The colonization and adaptation into new habitat can be 
observed in the star-like shape of network haplotype. 

Conclusions 

In conclusion, Lu. gomezi is a species with a higher genetic 
pool or variability and mild population structure, due to 
possible capacity migration and local adaptation to envir- 
onmental changes or colonization potential. However the 
existence of geographic barriers, such as the central 
mountain range separates a subpopulation. The establish- 
ment of Lu. gomezi in the Isthmus of Panama could be 
interpreted by the exchange of mammals in the American 
continent, vegetation and climate conditions of the epoch. 
Their relevance as a carrier of cutaneous leishmaniasis in 
Panama is important for health schemes, especially if we 
consider that sandflies have a wide-spread geographic 
distribution in Panama, and this distribution is possibly 
associated to the landscape changes caused by deforest- 
ation [20]. In this study, mitochondrial and nuclear analysis 
has provided important information to allow assembles 
about genetic population and evolutionary history useful to 
understand the implications of different population 
genetic structures for cutaneous leishmaniasis epidemi- 
ology. Thereby, proposing new perspectives for the 
control of leishmaniasis vectors. 

Additional files 
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Abbreviations 

CL: Cutaneous leishmaniasis; CB3-NIN: Cytb gene, inter-genetic region 
IGS-1; tRNA-Ser: Inter-genetic region IGS-2, and start of the NADH1; EF 
a-1: Elongation factor alpha-1. 

Competing interests 

The authors declare that they have no competing interests. 
Authors' contributions 

Conceived and designed the experiments: AVC, MGT. Performed the 
experiments: AVC. Analyzed the data: AVC. Contributed reagents/materials/ 
analysis tools/specimens identifications: MGT, JDAF. Wrote the paper: AVC, 
JDAF and MGT. All authors read and approved the final version of the 
manuscript. 

Acknowledgments 

We are grateful to Maykon Cristiano Passos, for helping with laboratory and 
analyses data; Jorge Edergarm Santos, Gustavo Martins and Juliana Lopes 
Fietto for comments in previous drafts of this manuscript; furthermore we 
thank the anonymous reviewers and editor for their comments. 

Funding 

This research was supported by the Secretaria Nacional de Ciencias y 
Tecnologia Innovadora (SENACYT) Grant Col09-008. 

Author details 

department of Medical Entomology, Instituto Conmemorativo Gorgas de 
Estudios de la Salud, Panama, Panama, department of Biology, Universidade 
Federal de Vigosa, Vigosa, MG, Brasil. 3 Centro de Referenda Nacional e 
Internacional para Flebotomineos/Colegao de Flebotomineos, Instituto Rene 
Rachou-Fiocruz, Belo Horizonte, MG, Brasil. 

Received: 12 July 2013 Accepted: 3 January 2014 
Published: 8 January 2014 

References 

1. Christensen HA, Herrer A, Telford SR: Enzootic cutaneous leishmaniasis in 
eastern Panama. II. Entomological investigations. Ann Trop Med Parasitol 
1972, 66:55-66. 

2. Christensen H, Fairchild GB, Herrer A, Johnson C, Young D, de Vasquez A: 
The ecology of cutaneous, leishmaniasis in the Republic of Panama. 

J Med Entomol 1 983, 20:463-484. 

3. Christensen H, Johnson C, Vasquez AM: Leishmaniasis cutanea en Panama: 
Un breve resumen. Rev Med Panama 1984, 9:182-1987. 

4. Darling ST: Oriental sore in Panama. Proc Canal Zone Med Assn 1910, 
2:7-20. 

5. Miranda A, Carrasco R, Paz H, Pascale JM, Samudio F, Saldaha A, Santamaria G, 
Mendoza Y, Calzada JE: Molecular epidemiology of American tegumentary 
leishmaniasis in Panama. Am J Trop Med Hyg 2009, 81:565-571. 

6. Herrer A, Telford S, Christensen HA: Enzootic cutaneous leishmaniasis in 
eastern Panama I: investigation of the infection among forest mammals. 
Ann Trop Med Parasitol 1971, 65:349-358. 

7. Tesh RB, Chaniotis BN, Aronson MD, Johnson KM: Natural host preferences 
of Panamanian Phlebotominae sandflies as determined by precipitin 
test. Am J Trop Med Hyg 1971, 20:150-156. 

8. Herrer A, Christensen HA: Epidemiological patterns of cutaneous 
leishmaniasis in Panama. Am J Trop Med Hyg 1973, 25:54-58. 

9. Saldana A, Chaves LF, Rigg CA, Wald C, Smucker JE, Calzada JE: Clinical 
cutaneous leishmaniasis rates are associated with household Lutzomyia 
gomezi, Lu. Panamensis, and Lu. trapidoi Abundance in Trinidad de Las 
Minas, Western Panama. Am J Trop Med Hyg 2013, 88(3):579. 

10. Jaramillo C, Travi BL, Montoya J: Vector competence of some neotropical 
sandflies for the Leishmania {Viannia) braziliensis complex. Med Vet 
Entomol 1994, 8:1-7. 

11. Jimenez AE, Rojas JC, Vargas F, Herrero MV: Temporal and spatial variation 
of phlebotomine (Diptera: Psychodidae) community diversity in a 
cutaneous leishmaniasis endemic area of Costa Rica. J Med Entomol 2000, 
37:216-221. 

12. Bejarano EE, Uribe S, Rojas W, Velez ID: Phlebotomine sand flies 
(Diptera: Psychodidae) associated with the appearance of urban 



Valderrama et al. Parasites & Vectors 2014, 7:9 
http://www.parasitesandvectors.eom/content/7/1/9 



Page 10 of 11 



leishmaniasis in the city of Sincelejo, Colombia. Mem Inst Oswoldo Cruz 
2002, 97:645-647. 

13. Santamana E, Ponce N, Zipa Y, Ferro C: Presencia en el peridomestic de 
vectores infectados con Leishmania {Viannia) panamensis en dos focos 
endemicos en el occidente de Boyaca, Piedemonte del valle del 
Magdalena medio, Colombia. Biomedica 2006, 26:82-94. 

14. Vivero RJ, Contreras MA, Cadena H, Acosta LA, Mondragon K, Velez A, 
Uribe SI, Velez ID: Sand flies Phlebotominae (Diptera: Psychodidae) in five 
Central American countries. In International Symposium on Phlebotominae 
Sandflies ISOPS VII: 25-30 april. Kudasy-Turkey; 201 1:95. 

15. Martins AV, Williams P, Falcao AL: American sand flies (Diptera: Psychodidae: 
Phlebotominae). Rio de Janeiro: Acad Bras Gene; 1978:195. 

16. Young DG, Duncan MA: Guide to the identification and geographic 
distribution of Lutzomyia sandflies in Mexico, the West Indies, Central 
and South America (Diptera: Psychodidae). Mem Amer Ent Inst 1994, 
54:1-881. 

17. Galati EAB: Classificacao de Phlebotominae. In Flebotomfneos do Brasil. 
Edited by Rangel ER, Lainson R. Rio de Janeiro: Editora Fiocruz Brazil; 
2003:23-52. 

18. Chaniotis BN, Correa MA, Tesh RB, Jonhson KM: Daily and seasonal 
man-biting activity of Phlebotominae sandflies in Panama. J Med Entomol 
1971,8:415-420. 

19. Travi BL, Alder GH, Lozano M, Cadena H, Montoya-Lerma J: Impact of 
habitat degradation on Phlebotominae (Diptera: Psychodidae) of tropical 
dry forest in northern Colombia. J Med Entomol 2002, 39:451-456. 

20. Valderrama A, Tavares MG, Andrade Filho DJ: Anthropogenic influence on 
the distribution, abundance and diversity of sandfly species (Diptera: 
Phlebotominae: Psychodidae), vectors of cutaneous leishmaniasis in 
Panama. Mem Inst Oswaldo Cruz 201 1, 106:1024-1031. 

21. Feliciangeli MD, Rabinovich J: Abundance of Lutzomyia ovallesi but not 
Lu. gomezi (Diptera: Psychodidae) correlated with cutaneous 
leishmaniasis incidence in north-central Venezuela. Med Vet Entomol 
1988, 12:121-131. 

22. Duque P, Velez ID, Morales M, Sierra D: Sand flies fauna involved in the 
transmission of cutaneous leishmaniasis in Afro-Colombian and 
Amerindian communities of Choco, Pacific Coast of Colombia. 
Neotrop Entomol 2004, 33:255-264. 

23. Cortes LA, Fernandez JJ: Especies de Lutzomyia en un foco urbano de 
leishmaniasis visceral y cutanea en El Carmen de Bolivar, Bolivar, 
Colombia. Biomedica 2008, 28:433-440. 

24. Jorquera A, Gonzalez R, Marchan-Marcano E, Oviedo M, Matos M: Multiplex- 
PCR for detection of natural Leishmania infection in Lutzomyia spp 
captured in an endemic region for cutaneous leishmaniasis in state of 
Sucre, Venezuela. Mem Inst Oswaldo Cruz 2005, 1 00:43-46. 

25. Vivero RJ, Contreras-Gutierrez MA, Bejarano EE: Analisis de la estructura 
primaria y secundaria del ARN de transferencia mitocondrial para serina 
en siete especies de Lutzomyia. Biomedica 2007, 27:429-438. 

26. Azpurua J, de la Cruz D, Valderrama A, Windsor D: Lutzomyia sandfly 
diversity and rates of infection by Wolbachia and an exotic Leishmania 
species on Barro Colorado Island, Panama. PLoS Negl Trop Dis 2010, 4:e627. 

27. Ready PD, Day JC, de Souza AA, Rangel EF, Davies CR: Mitochondrial DNA 
characterization of populations of Lutzomyia whitmani (Diptera: 
Psychodidae) incriminated in the peridomestic and silvatic transmission 
of Leishmania species. Bull Entomol Res 1 997, 87:1 87-1 95. 

28. Soto SI, Lehmann T, Rowton ED, Velez BID, Porter CH: Speciation and 
population structure in the morphospecies Lutzomyia longipalpis (Lutz 
& Neiva) as derived from the mitochondrial ND4 gene. Mol Phylogen Evol 
2001, 18:84-93. 

29. Hodgkinson VH, Birungi J, Quintana M, Deitze R, Munstermann LE: 
Mitochondrial Cytocrome B variation in population of the visceral 
leishmanisis vector Lutzomyia longipalpis across eastern Brazil. Am J Trop 
Med Hyg 2003, 69:386-392. 

30. Hamarsheh O, Presber W, Abdeen Z, Sawalha S, Al-Lahem A, Schoenian G: 
Genetic structure of mediterranean populations of the sandfly 
Phlebotomus papatasi by mitochondrial cytochrome b haplotype 
analysis. Med Vet Entomol 2007, 21:270-277. 

31. Esseghir S, Ready PD, Killick-Kendrick R, Ben-Ismail R: Mitochondrial 
haplotypes and phylogeography of Phlebotomus vectors of Leishmania 
major. Insect Mol Biol 1 997, 6:21 1 -225. 

32. Depaquit J, Lienard E, Verzeaux-Griffon A, Ferte A, Bounamous H, Gantier JC, 
Hanafi HA, Jacobson RL, Maroli M, Moin-Vaziri V, Muller F, Ozbel Y, 



Svobodova M, Volf P, Leger N: Molecular homogeneity in diverse 
geographical populations of Phlebotomus papatasi (Diptera, Psychodidae) 
inferred from ND4 mtDNA and ITS2 rDNA Epidemiological consequences. 

Infect Genet Evol 2008, 8:1 59-1 70. 

33. Zeledon R, Hidalgo H, Viquez A, Urbina A: Atypical cutaneous 
leishmaniasis in a semiarid region of north-west Costa Rica. Trans R Soc 
Trop Med Hyg 1989, 83:786. 

34. Warburg A, Saraiva E, Lanzaro GC, Titus R, Neva F: Saliva of Lutzomyia 
longipalpis sibling species differs in its composition and propensity to 
enhance leishmaniasis. Phil Trans Roy Soc Lond 1994, 345:223-230. 

35. Lourenco-de-Oliveira R, Vazeille M, Fillipis AMB, Failloux A: Aedes aegypti in 
Brazil: genetically differentiated populations with high susceptibility to 
dengue and yellow fever viruses. Tran R Soc Trop Med Hyg 2004, 98:43-54. 

36. Black WC IV, Tabachnick WJ: Population genetics in vector biology. In The 
Biology of Disease Vectors. Edited by Beaty BJ, Marquardt WC. Colorado: 
University Press of Colorado, Niwot; 2004:417-437. 

37. Reisen WK, Lothrop HD, Presser SB, Hardy JL, Gordon EW: Landscape 
ecology of arboviruses in southeastern California: temporal and spatial 
patterns of enzootic activity in Imperial Valley, 1991-94. J Med Entomol 
1997, 34:179-188. 

38. Sudia WD, Chamberlain RW: Battery operated light trap, an improved 
model. Mosq News 1962, 22:126-129. 

39. Michalsky EM, Fortes-Dias CL, Pimenta PFP, Secundino NFC, Dias ES: 
Assessment of PCR in the detection of Leishmania spp in experimentally 
infected individual phlebotomine sandflies (Diptera: Psychodidae: 
Plebotominae). Rev Inst Med Trop Sao Paulo 2002, 44:255-259. 

40. Testa JM, Montoya-Lerna J, Cadena H, Oviedo M, Ready PD: Molecular 
identification of vectors of Leishmania in Colombia: mitochondrial 
introgression in the Lutzomyia townsendi series. Acta Trop 2002, 
84:205-218. 

41. Parvizi P, Assmar M: Nuclear elongation factor 1-a gene, A molecular 
marker for Iranian sandfly identification. Iranian J Publ Health 2007, 
36:25-37. 

42. Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: 
Molecular evolutionary genetics analysis using maximum likelihood, 
evolutionary distance, and maximum parsimony methods. Mol Biol Evol 
2011,28:2731-2739. 

43. Librado P, Rozas J: DnaSP v5: A software for comprehensive analysis of 
DNA polymorphism data. Bioinformatics 2009, 25:1451-1452. 

44. Excoffier L, Laval G, Schneider S: Arlequin ver. 3.0: An integrated software 
package for population genetics data analysis. Evol Bioinform Online 2005, 
1:47-50. 

45. Weir BS, Cockerham C: Estimating F-statistics for the analysis of 
population structure. Evolution 1984, 38:1358-1370. 

46. Excoffier L, Smouse PE, Quattro JM: Analysis of molecular variance inferred 
from metric distances among DNA haplotypes: application to human 
mitochondrial DNA restriction data. Genetics 1992, 131:479-491. 

47. Bandelt HJ, Forster P, Rohl A: Median-joining networks for inferring 
intraspecific phylogenies. Mol Biol Evol 1999, 16:37-48. 

48. Huelsenbecr JP, Ronquist, MrBayes F: Bayesian inference of phylogenetics 
tree. Bioinformatics 2001, 17:754-755. 

49. Nylander JAA, Ronquist F, Huelsenbeck JP, Nieves-Aldrey JL: Bayesian 
phylogenetic analysis of combined data. System Biol 2004, 53:47-67. 

50. Tajima F: Evolutionary relationship of DNA sequences in finite 
populations. Genetics 1983, 105:437-460. 

51 . Fu YX: Statistical tests of neutrality of mutations against population 
growth, hitchhiking and background selection. Genetics 1997, 
147:915-925. 

52. Rogers AR, Harpending HC: Population growth makes waves in the 
distribution of pairwise genetic differences. Mol Biol Evol 1992, 
9:552-569. 

53. Brower AVZ: Phylogeny of Heliconius butterflies inferred from 
mitochondrial DNA sequences (Lepidoptera: Nymphalidae). Mol 

Phylogenet Evol 1 994, 3:1 59-1 74. 

54. Galtier N, Nabholz B, Glemin S, Hurst GD: Mitochondrial DNA as a marker 
of molecular diversity: a reappraisal. Mol Ecol 2009, 18:4541-4550. 

55. Bejarano EE, Rojas W, Uribe S, Velez ID: Genetic analysis of a recently 
detected urban population of Lutzomyia evansi (Diptera: Psychodidae) in 
Colombia. Rev Soc Entomol Argent 2009, 68:1 35-1 41 . 

56. Hedrick PW: Sex differences in mutation, recombination, selection, gene 
flow, and genetic drift. Evolution 2007, 61(12):2750-2771. 



Valderrama et al. Parasites & Vectors 2014, 7:9 
http://www.parasitesandvectors.eom/content/7/1/9 



Page 11 of 1 1 



57. Calcagnotto D: Taxa de evolucao e o relogio molecular. In Biologia 
Molecular e Evolucao. Edited by Matioli SR. Ribeirao Preto Sao Paulo: Hobs 
Editora; 2001:52-63. 

58. Nahum LA: Evolucao dos genomas. In Biologia Molecular e Evolucao. Edited 
by Matioli SR. Ribeirao Preto SP: Holos Editora; 201 1:82-96. 

59. Simon C, Frati F, Beckenbach A, Crespi B, Liu H, Flook P: Evolution, 
weighting and phylogenetic utility of mitochondrial gene sequences 
and a compilation of conserved polymerase chain reaction primers. 
Ann EntomolSocAm 1994, 87:651-701. 

60. Alexander JB, Young DG: Dispersal of phlebotomine sand flies (Diptera: 
Psychodidae) in a Colombian focus of Leishmania {Viannia) braziliensis. 
Mem Inst Oswaldo Cruz 1 992, 87:397-403. 

61 . Morrison AC, Ferro C, Morales A, Tesh RB, Wilson ML: Dispersal of the sand 
fly Lutzomyia longipalpis (Diptera: Psychodidae) at an endemic focus of 
visceral leishmaniasis in Colombia. J Med Entomol 1993, 30:427-435. 

62. Marquez LM, Lampo M, Rinaldi M, Lau P: Gene flow natural and domestic 
populations of Lutzomyia longipalpis (Diptera: Psychodidae) in a 
restricted focus of American visceral leishmaniasis in Venezuela. 

J Med Entomol 2001,38:12-16. 

63. De Queiroz BV, Coutinho-Abreu IV, Sonoda IV, Melo MA, de Andrade PP, 
de Castro JA, Rebelo JM, Carvalho SM, Ramalho-Ortigao M: Genetic 
structure of natural populations of the sand fly Lutzomyia longipalpis 
(Diptera: Psychodidae) from the Brazilian northeastern region. Acta Trop 
2006, 98(1 ):1 5-24. 

64. Rocha LS, Falqueto A, dos Santos CB, Ferreira AL, da Graea CC, Grimaldi G, 
Cupolillo E: Survey of natural infection by Leishmania in sand fly species 
collected in southeastern Brazil. Trans R Soc Trop Med Hyg 2010, 
104:461-466. 

65. Anez N, Cazorla D, Nieves E, Chataing B, Castro M, de Yarbuh AL: 
Epidemiologia de la leishmaniasis tegumentaria en Merida, Venezuela. I. 
Diversidad y dispersion de especies flebotominas en tres pisos 
altitudinales y su posible rol en la transmision de la enfermedad. 

Mem Inst Oswaldo Cruz 1 988, 83:455-463. 

66. Molina-Cruz A, de Merida AM, Mills K, Rodriguez F, Schoua C, Yurrita MM, 
Molina E, Palmieri M, Black WC IV: Gene Flow Among Anopheles albimanus 
populations in Central America, South America and the Caribbean 
Assessed by microsatellites and mitochondrial DNA. Am J Trop Med Hyg 
2004,71:350-359. 

67. Loaiza JR, Bermingham E, Scott ME, Rovira JR, Conn JE: Species 
composition and distribution of adult Anopheles (Diptera: Culicidae) in 
Panama. J Med Entomol 2008, 45:841-851. 

68. ETESA. http://www.hidromet.com.pa/. 

69. Condit R, Hubbell SP, Foster RB: Changes in a tropical forest with a 
shifting climate: results from a 50 ha permanent census plot in Panama. 

JTropEcol 1996, 12:231-256. 

70. Ready PD, de Souza AA, Macario Rebelo JM, Day JC, Silveira FT, 
Campbell-Lendrum D, Davies CR, Costa JML: Phylogenetic species and 
domesticity of Lutzomyia whitman! at the southeast boundary of 
Amazonian Brazil. Trans R Soc Trop Med Hyg 1998, 92:159-160. 

71 . Coates AG, Obando JA: Geological evolution of the Central American 
Isthmus. In Evolution and environment in tropical America. Edited by 
Jackson JBC, Budd AF, Coates AG. Chicago: University of Chicago Press; 
1996:21-56. 

72. Cox CB: Plate tectonics, seaways and climate in the historical 
biogeography of mammals. Mem Inst Oswaldo Cruz 2000, 95:509-516. 

73. Webb SD: The great American biotic interchange: patterns and 
processes. Ann Mo Bot Gard 2006, 93:245-257. 

74. Pehalver E, Grimaldi D: Assemblages of mammalian hair and 
blood-feeding midges (Insecta: Diptera: Psychodidae: Phlebotominae) in 
Miocene amber. Trans R Soc Edin: Earth Sci 2005, 96:1 77-195. 

75. Heckadon MS: Panama: Puente biologico. Panama. Panama: Smithsonian 
Tropical Research Institute; 2001:233. 

76. Hoernle KA, Bogaard PVD, Werner R, Lissinna B, Hauff G: Missing history 
(16-71 Ma) of the Galapagos hotspot: Implications for the tectonic and 
biological evolution of the Americas. Geology 2002, 30:795-798. 

77. Linares OF, Sheets PD, Rosenthal EJ: Prehistoric agriculture in tropical 
highlands. Science 1975, 187:137-145. 



78. 



79. 



Piperno DR, Jones J: Paleoecological and archaeological implications of a 
late pleistocene/early holocene record of vegetation and climate from 
the Pacific coastal plain of Panama. Quaternary Res 2003, 59:79-87. 
Marshall LG: Land Mammals and the Great American Interchange. Am Sci 

1988, 76(4):380-388. 

Ravelo AC, Andreasen DH, Lyle M, Lyle AO, Wara MW: Regional climate 
shifts caused by gradual global cooling in the Pliocene epoch. 

Nature 2004, 429:263-267. 



doi:1 0.1 186/1 756-3305-7-9 

Cite this article as: Valderrama et al.: Phylogeography of the Lutzomyia 
gomezi (Diptera: Phlebotominae) on the Panama Isthmus. Parasites & 
Vectors 2014 7:9. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at 
www.biomedcentral.com/submit 



o 



BioMed Central 



